Non-equilibrium transport through a point contact in the v — 5/2 non-Abelian quantum Hall state 



Adrian Feiguin, 1 ' 2 Paul Fendley, 3 Matthew P.A. Fisher, 2,4 and Chetan Nayak 2 ' 4 

1 Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park, MD 20742, USA 
Microsoft Research, Station Q, CNSI Building, University of California, Santa Barbara, California 93106, USA 
3 Department of Physics, University of Virginia, Charlottesville, VA 22904, USA 
4 Department of Physics, University of California, Santa Barbara, California, 93106, USA 

We analyze charge-e/4 quasiparticle tunneling between the edges of a point contact in a non-Abelian model 
of the v = 5/2 quantum Hall state in the presence of a finite voltage difference using the time-dependent density- 
matrix renormalization group (DMRG) method. We confirm that, as the voltage decreases, the system is broken 
into two pieces. In the limits of small and large voltage, we recover the results expected from perturbation theory 
about the infrared and ultraviolet fixed points. We test our methods by finding the analogous non-equilibrium 
current through a point contact at v = 1/3. 



Introduction- The v = 5/2 fractional quantum Hall state 
|D1 S is suspected to be a non-Abelian topological state. 
While theoretical evidence has been steadily accumulating 
over the years, there has been little experimental evidence - 
until now. Several transport measurements utilizing a point 
contact, one a low-frequency noise ('shot noise') measure- 
ment yj, the other the tunneling current [5], indicate that the 
smallest quasiparticle charge at this plateau is e/4. This is 
consistent with two non-Abelian models of the v = 5/2 frac- 
tional quantum Hall state JfJ 0, |H 01, but it also consistent 
with an Abelian model 

One limit to the success of these measurements is that the 
data is compared only to lowest-order perturbative calcula- 
tions valid for small inter-edge tunneling at the point contact. 
However, inter-edge tunneling in these experiments is not so 
small, and may have effects not described with low-order per- 
turbation theory. Therefore, it is important to compute the 
expected I — V curve (and the low-frequency noise) beyond 
the perturbative regime. In so doing, we follow the crossover 
from the limit of weak quasiparticle tunneling across the point 
contact to the low-temperature, low-voltage regime. 

In this paper, we numerically compute the zero-temperature 
current through a point contact in the presence of a finite volt- 
age bias in two fractional quantum Hall states. The first is the 
Abelian Laughlin state at v = 1/3, which allows us to check 
our numerical procedure, and to confirm that indeed the Bethe 
ansatz computations [11] are applicable out of equilibrium. 
The second is the non-Abelian Moore-Read Pfaffian state with 
v = 5/2 We show that in both cases the droplet eventu- 
ally breaks in two at low voltages. In the weak quasiparticle 
tunneling limit, the current follows the predicted power law. 
In the low-voltage limit, the conductance is obtained and is 
approached in a manner consistent with predictions lfl2ll . 

We use the time-dependent Density Matrix Renormaliza- 
tion Group (DMRG) method 10 Q, because other ap- 
proaches have difficulties. Namely, at finite bias, this non- 
equilibrium calculation is not amenable to a Monte Carlo sim- 
ulation. (The Monte Carlo computation in the v = 1/3 case 
is applicable only to the linear-response regime Jl5ll .) As we 
discuss below, the problem can be mapped onto resonant tun- 
neling between attractive Luttinger liquids. Since the 'leads' 



are interacting, a Wilsonian numerical renormalization group 
cannot be used, unlike in the Kondo problem. Finally, when 
v = 5/2, the model is not integrable, and thus the Bethe 
ansatz is not applicable. 

Models- The edge excitations of the Laughlin state at 
v = 1/3 are described by a chiral Luttinger liquid. At a point 
contact, two edges come into close proximity so that a charge 
e /3 quasiparticle can be backscattered from one to the other. 
We describe each edge by a chiral boson <f>i, so that a charge 
e/3 Laughlin quasiparticle is created by the operator e 1 *^/^. 
The effective Hamiltonian is 

H i/ 3= y f dx{dx(Pi) 2 + te i(*i(o)-^(o»/VS + hx 

i=1.2 J 

(1) 

The tunneling amplitude has scaling dimension [t] = 2/3. 
Hence, at zero temperature the backscattere current Ib ~ 
t 2 V" -1 / 3 in the limit of small tunneling current and large volt- 
ages. In the low voltage limit, perfect backscattering occurs 
and the Hall bar effectively breaks in two Charge trans- 
port between the two halves is due to electron tunneling, so for 
small V, Ib — ^ \ V ~ V . The picture was confirmed by 
finding the full crossover from weak to strong backscattering 
via a Monte Carlo calculation of the linear-response current at 
non-zero temperature [1511 and the Bethe ansatz solution for 
the full I B (V) curve lUl 

To find the zero-temperature Ib(V) curve numerically, we 
rewrite this problem as non-resonant tunneling between two 
semi-infinite non-chiral spinless Luttinger liquids. We define 
4> a (x) on the half-line x < and qbb (x) on the half-line x > 
as follows: 4> a {x < 0) = 4>i(x) + qt>i(—x) and qt>b{x > 0) = 
4>2(x) + 4>2{—x). The tunneling term in ((TJ then becomes 

H non — res — 

t (Vi(O)^(Q) +h.c) (2) 

where ip a t,(x) — e 4 *"' 6 ^^ 2 ^ are the two Luttinger quasipar- 
ticle creation operators (which are non-local combinations of 
the Laughlin quasiparticles). Since these operators have scal- 
ing dimension 1/6, the Luttinger liquids have g = 3 in the 
conventions of Ref. Ila a duality transformation maps this to 
a g = 1/3 Luttinger liquid perturbed by a <5-function impurity. 
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In our DMRG computations, we use a tight-binding model 
of spinless fermions with nearest-neighbor attractive interac- 
tions to describe each Luttinger liquid. This model is equiv- 
alent, under a Jordan- Wigner transformation, to an XXZ spin 
chain. Coupling the two liquids corresponds to including a 
link between the sites at the two ends, as illustrated in fig- 
ure la. Just as charge tunneling violates charge conservation 
of the individual edges, coupling the two chains violates the 
conservation of the individual magnetizations. We match pa- 
rameters by noting that the scaling dimension of the staggered 
spin-raising operator in the spin chain, Sl(—l) x , is equal to 
the Luttinger parameter g. The ferromagnetic XXZ spin chain 
anisotropy is then related to g by J z / J± = — cos(7r jig) lfl7ll . 
Thus iJ non _ rcs can be realized with J z /J±_ = — V3/2. 

Tunneling through a point contact for v = 5/2 also can be 
realized via coupled XXZ chains, by following the bosoniza- 
tion procedure of 11211 . We assume that the N = Landau 
level (of both spins) is filled and the N = 1 Landau level is 
in the half-filled Moore-Read Pfaffian state. The former are 
integer quantum Hall edge modes, and are the outermost ex- 
citations of the system; we ignore them because we focus on 
tunneling across the interior of a Hall droplet. The gapless 
chiral theory describing the edge excitations of the Moore- 
Read state consists of a free boson (the charge sector) and a 
free Majorana fermion (the neutral sector). We study the inter- 
edge backscattering of the basic charge e/4 quasiparticle at a 
point contact at x = 0. Charge-e/2 quasiparticles can also 
tunnel, as can neutral quasiparticles, but the latter does not 
affect the electrical conductivity and the former is expected 
to be smaller. The effects of the latter, as well as the exten- 
sion to the anti-Pfaffian and (3, 3, 1) states, will be discussed 
elsewhere 111 811 . 

To bosonize this model requires a fairly elaborate compu- 
tation, because the charge-e/4 quasiparticle has non-abelian 
statistics. When the dust settles, the tunneling Hamiltonian 
can be written in terms of two bosons <\> p and </> CT , and a Kondo 



spin S. The resulting Hamiltonian is 112 ] 



n 



5/2 



dx 



\7T i d x<t>p) 
\ Z7T 

+ t (Ve-^(°)/ 2 + S~e^ (0)/2 ) cos(0 p (O)/2) . (3) 



The tunneling amplitude has scaling dimension [t] = 3/4. 
Therefore, in the limit of small current, Ib ~ t 2 V~ x l 2 at 
T = 0; this limit occurs for large voltage. In the oppo- 
site limit, V — ► 0, inter-edge tunneling becomes strong, and 

i 2 

Ib — * g \ V. Deviations from this total backscattering vary 

i 2 c 

generically as: Ib — 5 \ V ~ V . However, when the tun- 
neling of neutral Majorana fermions is neglected, the system 
flows to the infrared fixed point along a special direction, so 
that/B — 4 t~ V ~ V 15 11211 . (Even though Majorana fermion 
tunneling does not directly contribute to charge transport, it 
affects the flow into the infrared fixed point.) In this paper, we 
will compute Ib(V) for arbitrary V. 

The Hamiltonian (01 has the form of resonant tunneling 



between attractive Luttinger liquids; the reason that the Lut- 
tinger liquids are attractive (g > 1) is that the tunneling oper- 
ator has scaling dimension 1/4, which is highly relevant. We 
can make the relation to resonant tunneling between Luttinger 
liquids more apparent by rewriting (f3]l in the form 



dx 



T~ \ {d x (j)a) 2 + {d x (j> b y 

zir \ 



+ t (ft e **.(o)A/5 + t it e <*»(o)A/ff + h .c. (4) 



where S = d, S + = d' annihilate/create a particle on the 
resonant level. The Luttinger coupling is g = 2 with <\> a jb = 

As with the earlier case, we utilize two semi-infinite ferro- 
magnetic XXZ chains for the two Luttinger liquids. Here, the 
two chains couple through the resonant level, which in the lat- 
tice model corresponds to adding an extra site, as illustrated 
in figure lb. Since g = 2 here, we have J z /J± = — l/v2- 

We apply a potential difference V between the two leads, 
or potential differences ±V/2 between lead a/b and the 
resonant level. This modifies the tunneling term in (O to 
t rf e i<f>a/V5 e ieVt/2 + td \ e ici >b /^g e -teVt/2 + hc However, 

in the original 5/2 point contact problem, the tunneling term 
transfers charge e/4 between the two edges, so a potential 
difference V between the edges modifies the tunneling term 
in ® to t (S+e-^/ 2 + S^e' 1 ^/ 2 ) cos(0 p (O)/2 + eVt/4). 
Because the current is proportional to the charge squared, the 
relation between the current in the original 5/2 point contact 
problem and in the Luttinger liquid resonant tunneling prob- 
lem is /b MR = (1/2) 2 ires. Likewise in the non-resonant case 
appropriate for v = 1/3, we have Ib 1 ^ = (l/3) 2 / n on-rcs- 

Time-Dependent DMRG- We find the I-V curves of the 
non-resonant and resonant tunneling problems using the time- 
dependent DMRG method. Two fundamental aspects of our 
calculation make it particularly unsuited to more conventional 
techniques such as Wilson's numerical renormalization group: 

(i) the conducting leads are interacting Luttinger liquids and 

(ii) we are interested in the non-linear regime, i.e. large volt- 
age bias. In ID metallic systems, correlations can drastically 
affect the density of states and transport properties. In par- 
ticular, repulsive interactions suppress charge transport, while 
attraction "heals" the system, enhancing the conductance. The 
time-dependent DMRG is well suited for our systems, be- 
cause it allows one to seamlessly incorporate interactions into 
the leads, and is not restricted to the linear response regime. 

Our technique consists of evaluating the time dependence 
of the current though the weak link or quantum dot, after a 
voltage bias is applied [19, 23]. In a first step, the ground state 
is calculated using the conventional DMRG technique. Then, 
the system is quenched: by applying a shift in the chemical 
potential S/j,l = V/2, 5^lr = —V/2 to the left and right 
leads, respectively. The resulting non-equilibrium system is 
evolved in time by solving the time-dependent Schrodinger 
equation. As a response to the quench, a current starts flowing 
through the system. Typically, the current grows and a tran- 
sient is observed in the beginning, followed by oscillations 
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FIG. 1: (a) Junction model used to study the current though a quan- 
tum point contact at v = 1 /3, consisting of a weak link with connect- 
ing two interacting leads, (b) Quantum dot, or resonant level, system 
used for v — 5/2. (c) Current through a junction after a step bias is 
applied. We show results for L = 120, t q = 0.1, V = 0.04, and dif- 
ferent boundary conditions, (d) Time-dependent current though the 
same system, with damped boundary conditions, and different values 
of the bias V . Time is measured in units of the hopping t q . 



that tends to stabilize at a constant value, corresponding to the 
steady state |20[ 21]. Since the leads used in the calculation 
are finite, a reversing of the current is observed after the wave 
packet reaches the boundaries and is reflected back. This de- 
termines a time scale in which we expect the current to sta- 
bilize at a plateau value. As we show here, depending on the 
choice of parameters in the model, it is sometimes necessary 
to study large systems in order to achieve a steady state. In 
some cases, when the transient region is large, and the system 
too small, this is hard to attain. In order to improve the behav- 
ior of the system, we used long leads and damped boundary 
conditions, by exponentially decreasing the coefficients in the 
Hamiltonian toward the end of the chains. As a result, remov- 
ing a particle from these regions becomes energetically costly, 
and they effectively behave as reservoirs. As a consequence, 
the charge becomes trapped and accumulates without getting 
reflected, leading to longer plateaus lfl9[|22ll . 

Results- The tunneling problem between i/ = 1/3 edges 
reduces to studying the current though a weak link, or junc- 
tion, connecting two interacting spinless leads, while the prob- 
lem at v = 5/2 corresponds to a resonant level, or quantum 
dot, as seen in FigsQIa) and (b), respectively. These two mod- 
els were studied in an early formulation of the time-dependent 
DMRG method [24]. The smaller the inter-lead hopping am- 
plitude t q , the larger is the initial transient in the current, mak- 
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FIG. 2: Ib-V characteristics of quasi-particle tunneling between 
v = 1/3 quantum Hall edges, as modelled by a junction system 
modeling, for different values of the tunneling amplitude t q , obtained 
using the time-dependent DMRG. Error bars correspond to the errors 
in averaging the current over an interval of time of the order of 40, in 
unit of the hopping, and system sizes up to L — 160. Lines are fits 
to the data using the exact Bethe ansatz solution. 



ing it more difficult to reach a steady state in a finite system. 
We found it necessary to use long chains, up to 160 sites. In 
FigtHc), we show results for v = 1/3, comparing the behav- 
ior of the current in systems with different sizes and boundary 
conditions. The damped boundary conditions, while yielding 
the same steady current for given values of t q and bias V, 
extend the duration of the plateau, allowing one to reach the 
steady state in smaller systems. In FigQId) we show results 
for t q = 0.1, and different values of the bias V, for a sys- 
tem with 120 sites and damped boundary conditions. Typical 
simulations extend to times of the order of 300 in units of the 
hopping, using a 3rd order Suzuki-Trotter decomposition of 
the evolution operator with a time step t = 0.2 and keeping 
the truncation error below 10~ 7 lr25ll . The current is averaged 
over an interval of time, and the error calculated following the 
prescription discussed in Refs. IJjl 23]. The Suzuki-Trotter 
error associated to the finite time-step was found to be much 
smaller than the error in the average. 

In Fig. |2] we show the Ib-V characteristic curves for the 
v = 1/3 case, for different values of t q . At small biases, the 
system exhibits a conductance G aon - Tea = 3e 2 /h which cor- 



responds to Gi/3 = |e /h, as expected from the arguments 
of Ref. Ill 611 . As the bias grows, the system departs from the 
linear response regime and crosses over to the scaling behav- 
ior associated with quasiparticle tunneling at the ultraviolet 
fixed point, I B oc V^ 1 / 3 ||I^] (note that we plot I B /V vs. 
V). As seen in Fig|2l not only the asymptotic power law but 
the full crossover follows the behavior predicted by the Bethe 
ansatz solution [11]. The exact Bethe ansatz expression has 
a free parameter, corresponding to the tunneling amplitude or 
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turn Hall droplet to split into two droplets which are coupled 
through weak electron hopping, as predicted in Ref.[12j. They 
further enable us to access the crossover regime of intermedi- 
ate biases where we find quantitative deviations from power- 
law behavior. Further work will compute IU8I1 the analogous 
I B -V curves for the anti-Pfaffian H H and (3, 3, 1) states 
ifioll and compare all three to experimental measurements JH. 
Such a comparison could pave the way to correctly identifying 
the v = 5/2 quantum Hall state. 
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FIG. 3: Ib-V characteristics of quasi-particle tunneling between 
v = 5/2 quantum Hall edges, as modelled by a quantum dot system, 
for different values of the tunneling amplitude t q , obtained using the 
time-dependent DMRG. Error bars correspond to the errors in aver- 
aging the current over an interval of time of the order of 40, in unit 
of the hopping, and system sizes up to L = 140. Lines are fits to the 
data using an expression Ib ~ V -1 ^ 2 . 

t q , which can be used to fit the numerical results. The agree- 
ment is excellent for small t q , up until the very large biases 
at which the lattice model no longer accurately represents the 
quantum Hall edge due to curvature of the dispersion. At large 
t q , the system remains near the infrared fixed point up to large 
biases, and the scaling behavior associated with quasiparticle 
tunneling at the ultraviolet fixed point cannot be observed. 

For the v = 5/2 case, we followed the same procedure 
described above, but using the resonant level system shown 
in FigQJb). The results are depicted in Fig. [3] At small 
bias, the system exhibits a conductance G Tes = 2e 2 /h which 
corresponds to Gmr = h^ 2 lh, as expected from the argu- 

| I " 

ments of Ref. 112J. As the bias is increased for small t q , the 
system crosses over to the asymptotic power-law associated 
with charge-e/4 quasiparticle tunneling, Ib oc V~ x / 2 . Again, 
at large bias, the numerical results depart from the universal 
regime, and exhibit the effects of the lattice. For small V, the 
deviations from G rcs = 2e 2 /h are too small to be reliably 
fit to a power law; this may be an indication that they are, in- 
deed ~ V 15 . In further work iflill . we will investigate whether 
the addition of an additional marginal operator will lead to the 
generic flow into this infrared fixed point with deviations from 
the asymptotic value ~ V 5 . 

Discussion. These results clearly demonstrate that weak 
inter-edge quasiparticle tunneling causes a Moore-Read quan- 
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